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ABSTRACT 

Using a sample of 89 snapshots from 58 hydrodynamic binary galaxy major merger simulations, we 
find that stellar remnants are mostly oblate while dark matter halos are mostly prolate or triaxial. 
The stellar minor axis and the halo major axis are almost always nearly perpendicular. This can 
be understood by considering the influence of angular momentum and dissipation during the merger. 
If binary mergers of spiral galaxies are responsible for the formation of elliptical galaxies or some 
subpopulation thereof, these galaxies can be expected to be oblate and inhabit their halos with the 
predicted shapes and orientations. These predictions are relevant to observational studies of weak 
gravitational lensing, where one must stack many optically aligned galaxies in order to determine the 
shape of the resulting stacked mass distribution. The simple relationship between the dark and lumi- 
nous matter presented here can be used to guide the stacking of galaxies to minimize the information 
lost. 

Subject headings: galaxies: formation — galaxies: interactions — galaxies: kinematics and dynamics 
— galaxies: structure 



1. INTRODUCTION 

The shapes and mass profiles of dark matter halos 
from cosmological iV-body simulations have long been 
studied (Dubinski fc CarlberglllQQlUNavarro et alJlloM 
lAllgQod et al.l l2006. and references therein). Cosmolog- 
ical simulations still lack sufficient resolution to track 
the shape and orientation of galaxies within their dark 
matter halos. There is no reason to believe that the 
shapes of galaxies and dark matter halos should be simi- 
lar. It has only recently become feasible to perform large 
suites of high-resoluti on bin ary galax y merger simula- 
tions (Naab & Burkcrt 2003; Cox 200l ICox et al.ll2?ffl5l: 
iRobertson et al..2006.) . and wc here use such simulations 
in order to study the shapes of the resulting galaxies and 
their host halos statistically. 

Observationally, the intrinsic shapes of elliptical galax- 
ies have remained elusive. It has long been known that 
there seem to be at least two classes of elliptical galaxies: 
massiv e, anisotropic galaxies ari d lower mass, oblate ro- 
tators fBendcr ct al."1988', '1992^. However, allowing the 
possibili ty of triaxia lity lea ds to degeneracie s in depro- 
iection iFranx et all 1X9911) . lAlam fc RvdenI t2002) and 
iVincent fc RvdenI l)2005f) have used Sloan Digital Sky 
Survey (SDSS) data to conclude that not all elliptical 
galaxies can be oblate. 

The relative orientations of galaxies and their dark 
halos is relevant to studies of weak gravitational lens- 
ing. Observers stack many images of galaxies in order to 
use the average deformation of the shapes of background 
galaxies to infer properties of the foreground mass dis- 
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tribution. It is important to stack galaxies coherently 
in order to build up a detectable signal. The model pre- 
sented here represents a physically well-motivated Ansatz 
to help interpret the results of weak lensing observations. 
Section[21gives a description of the galaxy merger simula- 
tions and our method of determining the shape of merger 
remnants, '^ gives our results, and SQl summarizes our 
conclusions. 

2. METHODS 

We analyze the shapes of 89 snapshots from 58 of these 
simulations. We study two samples of simulations. One 
is the "G" series, which consists of major and minor 
mergers with progenitor spiral galaxies typical of the 
nearby universe and spanning a factor of 40 in bary- 
onic mass and 20 total mass. In order to reduce the 
dependence on the progenitor galaxy model, here we only 
consider major mergers with mass ratios of 1:1 (G3-G3, 
G2-G2, Gl-Gl, and GO-GO) and roughly 3:1 (G3-G2, 
G2-G1, and Gl-GO). We also analyze the "Sbc" series of 
merger simulations, which are 1:1 major mergers of mas- 
sive, gas-rich spirals using a variety of different orbits 
and orientations. 

To calculate the shape of a merger remnant, we iter- 
atively diagonalize a moment of inertia tensor using an 
ellipsoidal window IjDubinski fc Carlberdll99l|) : 



My = Y^krrikn^krj^k 



(1) 



where r^^^ is the position vector, i,j refer to coordinates, 
and k re fers to parti c le num ber. The triaxial radius is 
given bv lFranx et alJ l|199lD : 

C = y/xya^+yyb^ + z^/c^ (2) 

where a, 6, and c are the major, intermediate, and minor 
axis lengths, respectively. The sum over k includes all 
particles for which r lies within the ellipsoid C — ^- The 
iteration is started with a spherical window {a = b = 
c — baryonic half-mass radius), and after each iteration 
a,b, and c are scaled so that half of the baryonic mass 
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Fig. 1. — Shapes of the luminous and dark components of simulated merger remnants. Left: the intermediate-to-major axis ratio vs. 
the minor-to-major axis ratio for stars. Right: Same as left panel, but for dark matter halos. Objects near the diagonal line are prolate 
spheroids, objects near b/a = 1 are oblate spheroids, and objects in between are triaxial. Dotted lines indicate constant triaxiality T. Most 
stellar remnants are oblate with e = 0.5. Most dark matter halos are either prolate or triaxial. 



is enclosed. The result does not appreciably change if 
equation (^ is modified to include C,'^ in the denominator. 
Using a spherical window rather than an ellipsoidal one 
results in systematically larger axial ratios but does not 
change the main result. 
Three-dimensional shapes of galaxies can be quantified 



with the triaxiality parameter T = {a — b )/{a 
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We call an object oblate, triaxial, or prolate if T < 0.25, 
0.25 < T < 0.75, or 0.75 < T, respectively. Shapes 
of galaxies can also be quantified by ellipticity e = 1 — 
b/a. EUipticities are most often used to describe two- 
dimensional shapes; we occasionally refer to the three- 
dimensional ellipticity of perfectly prolate or oblate {T = 
or 1) objects since there is no ambiguity about the use 
of the equation. 

Simulations were performed using the cntropy- 
conserving version of the S PH code GADGET 
l|Springel et all |200H iSorineel fc Hcrnauist 2002) with 
a gravitational smoothing length of 100 pc. The pro- 
genitor galaxies have baryonic masses from 1.6 x 10^ to 
2 X 10^"'^ M0, gas fractions between 20% and 70%, consist 
of ~100,000 particles, and use a parameterization of star 
formation feedback fro m supernovae tu n ed to match the 
empirical Schmid t law l|Kennicuttll998ri . iCoxl l|2004fl and 
ICox et aLl l|2005j) contain further information about the 
simulations. 
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Fig. 2. — The angle between the minor axis of the stars and 
the major axis of the dark matter. The stellar material is mostly 
oblate, the dark matter halos are mostly prolate, and the "pre- 
ferred" axes of the two shapes are nearly always perpendicular. 



3. RESULTS 

Figure ^ illustrates that most stellar remnants are 
oblate, while the dark matter halos in which they reside 
are mostly prolate or triaxial. Figure |2] shows that the 
short axis of the stellar system and the long axis of the 
dark matter halo are almost always nearly perpendicu- 
lar. This can be understood simply in terms of angular 
momentum and dissipation, as shown in Figure |31 

This model helps interpret the findings from studies of 
weak gravitational lensing. Hoekstra ct al. ( , 2004') find 
that the ellipticity of dark halos is 0.77lg'2i times the 
ellipticity of the light (i.e., halos are somewhat less flat- 
tened than galaxies), assuming that the two are aligned. 
According to our result, elliptical galaxies would either 
show an elliptical halo (if the long axis of the prolate halo 
is in the plane of the sky) or a circular halo (if the long 
axis of the halo is pointed toward the observer). Thus, 
the fiattening of the dark matter would follow that of the 
luminous matter, in agree ment with these o bservations. 

The interpretation of the IHoekstra et al.l data is com- 
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Fig. 3. — The physical interpretation of Figs. [Tlandl^in terms of 
angular momentum and dissipation. The total angular momentum 
in a merger simulation is usually dominated by the orbital angular 
momentum of the two galaxies. As the galaxies merge, both the 
luminous and dark components acquire angular momentum from 
the orbits of the progenitors. Their velocity dispersion increases 
along the axis parallel to the direction of approach, leading to an 
anisotropic velocity dispersion tensor. Gas in the simulation cools 
while largely conserving angular momentum so that it spins up to 
the point where the shape of the resulting stellar system is deter- 
mined by rotation, not velocity dispersion anisotropy. Meanwhile, 
the dark matter cannot cool, so its shape is determined by velocity 
dispersion anisotropy. Therefore, the stellar system is oblate, dark 
matter halos are prolate, and the angle between the "preferred" 
axes of these two shapes is ~90° . 
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mass density and po is a constant l|Navarro et al.llT99a: 
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Fig. 4. — Consequences of our model for weak gravitational 
lensing measurements where one must stack many optically aligned 
galaxies in order to detect flattening of the dark matter halo. The 
plot shows the apparent ellipticity of the halo mass surface den- 
sity vs. the apparent optical ellipticity of the stacked galaxies. 
Each shaded region shows the range of possible apparent ellip- 
ticities given an intrinsic, three-dimensional halo ellipticity. The 
lower bound of each shaded region is given by assuming a prolate 
halo and the upper bound by assuming an oblate halo. This is 
the result one would obtain if all galaxies followed the trend noted 
in this Letter and one stacked galaxies only with a given optical 
ellipticity. The observed halo ellipticity goes to zero for appar- 
ent optical ellipticities near zero because the galaxies cannot be 
oriented so that the stacking is coherent, even though each indi- 
vidual halo will have a nonzero projected ellipticity. The observed 
halo ellipticity only equals the three-dimensional ellipticity when 
halos are intrinsically oblate and the galaxies are viewed edge on; 
otherwise the flattening is underestimated. 



plicated by the inclusion of spiral as well a s elliptical 
galaxies in the sample. iMandelba um et al.l l)2006alhf) 
haye done a similar study using SDSS galaxies separated 
by Hubble type and found that the projected halo shapes 
for elliptical galaxies are aligned with the projected stel- 
lar shapes, in agreement with Hoekstra et al. Finally, the 
projected positions of satellite galaxies also seem to indi- 
cate that the proj ected shapes of elliptical galaxies and 
halos are aligned l|Sales fc Lambajl200l: lBrainerdll2005t 
lYang et alJl20nl . 

Weak lensing studies necessarily underestimate the 
flattening of dark matter halos. Figures01and|Slquantify 
this by simulating the weak gravitational lensing obser- 
vations. Given assumptions about the three-dimensional 
shapes and mass profiles of galaxies and their halos and 
a scheme for combining many galaxies into a single mass 
surface density, these two figures show shapes of the pro- 
jected halo mass surface densities. They allow observers 
to translate their two-dimensional measurements to a 
range of possibilities for the three-dimensional structure 
of dark matter halos. 

The hydrodynamic simulations discussed here do not 
represent a cosmologically unbiased sample, so they are 
not used as input to the simulated lensing observations. 
Instead we adopt a slightly idealized version of the cor- 
relation between halos and galaxies noted in this Letter. 
Nearly all of the baryonic components of the simulated 
galaxies are close to 2:1 oblate spheroids, so we assume 
that all early-type galaxies are so described. Thus, there 
is a simple mapping between viewing angle and optical 
ellipticity. We assume all galaxies follow the correla- 
tion between halos and galaxies noted here and that the 
halo mass density is given by a triaxial Navarro-Frenk- 
White profile: p = po/(C/^s)(l + C/fsY, where p is the 
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IContopoulosI 1)19561) showed that for a triaxial ellipsoid 
with constant three-dimensional axis ratios, the contours 
of constant projected surface density are ellipses with 
constant ellipticity and position angle, independent of 
the radial density profile. We only use the ellipticity and 
position angle of the baryonic component, so the radial 
profile of the baryons does not matter. TheEontopouloj 
()1956(1 analysis does not apply to the stacked dark matter 
halos, so Figures 0] and |5l depend on the radial density 
distribution of the halos. In practice the difference is not 
large. 

To simulate weak lensing measurements, we align the 
projected mass distributions based on projected light dis- 
tributions, stack the projected halo mass distributions, 
and fit an ellipse to the halo mass surface density dis- 
tribution where the area of the ellipse is constrained to 
equal 7r(3rs)^. This size for the ellipse is motivated by 
the approximate radius at which weak lensing observa- 
tions are sensitive to the halo shape (M. J. Hudson 2006, 
personal communication). The stacking either assumes 
a given inclination of the optical galaxy, averaging over 
the azimuthal angle (as in Fig. ^, or assumes that some 
minimum optical ellipticity is required to be included in 
the stack, averaging over the portion of the unit sphere 
that gives rise to sufficient optical ellipticities (as in Fig. 
0). 

Figure El shows that galaxies with low optical elliptic- 
ities will have low halo ellipticities because there is no 
preferred axis to use to stack galaxies. The only situa- 
tion where the projected halo ellipticity equals the three- 
dimensional halo ellipticity is when all stacked galaxies 
are viewed edge-on and halos are intrinsically oblate. 
Flattening is underestimated in all other cases. 

Figure shows the result of the more realistic scenario 
where all galaxies with optical ellipticities greater than 
some value are included in the stack. This allows one 
to transform projected ellipticities to three-dimensional 
ellipticities. For example, if an observer sets the min- 
imum optical ellipticity to 0.2 and measures a stacked 
halo ellipticity of 0.25, one can conclude that the three- 
dimensional ellipticity of halos is either 0.3 (for oblate 
halos), 0.5 (for prolate halos), or somewhere in between. 

As one enforces tighter constraints on the optical el- 
lipticity, the halo ellipticity goes up, but the cost is that 
fewer galaxies will make it into the stack. Under simple 
assumptions, one can estimate the signal-to-noise ratio 
(S/N) of the halo ellipticity measurement to be 



(S/N)tot = e2DV^^^^M (3) 

where €2Yy is the apparent ellipticity of the stacked halo 
mass surface density, $7 is the solid angle of viewing an- 
gles for which a galaxy will be included in the stack, N^q^ 
is the total number of galaxies in the survey, and txi is 
the error on the halo ellipticity when only one galaxy is 
used. We define ^ as the part of this expression which 
depends on the signal and the available solid angle: 



62DVf7 (4) 

Figure|Slthus also allows observers to estimate the quality 
of their measurement given the size of their survey and an 
estimate of the one-galaxy error on the halo ellipticity. In 
reality, observers do not know 17, but they could estimate 
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Fig. 5. — Weak gravitational lensing measurements when all galaxies with optical ellipticities greater than some value are included in 
the stacking. The shaded regions are defined as in Fig. |4] In both panels the x-axis is the minimum apparent optical ellipticity of galaxies 
included in the stacking. Left: The y-axis shows the apparent ellipticity of the stacked halo mass surface density. This is the halo ellipticity 
that would be measured if one stacked all galaxies with optical ellipticity greater than the given value. Right: Effect of different optical 
ellipticity requirements on the S/N ratio of the halo ellipticity measurement. Increasing the minimum optical ellipticity increases the signal 
(as seen in the right panel) but also increases the noise by reducing the number of galaxies included in the stack. The jz-axis is '3/ defined by 
eq. J3}. For small optical ellipticity cuts, the increase in signal is almost exactly canceled by the decrease in available solid angle, meaning 
that nothing is lost (or gained) by removing round galaxies from the stack. In reality, there is an advantage to removing low-ellipticity 
galaxies since their position angle is ill-defined. (This analysis assumes no error in position angle.) Above a minimum ellipticity of ~ 0.25 
the quality of the measurement declines rapidly because of smaller solid angle of available viewing directions. 



it from the from the minimum eUipticity of galaxies in 
their sample as long as our assumption that the stellar 
remnants are 2:1 oblate spheroids is not far wrong. 

4. DISCUSSION 

We have analyzed the three-dimensional shapes of 
galaxies and dark matter halos resulting from more than 
100 simulations of gas-rich galaxy mergers. Stellar rem- 
nants are nearly all oblate, with a few examples of triax- 
iality in the most gas-poor mergers. Dark matter halos 
are either prolate or triaxial, and the short axis of the 
baryons is perpendicular to the long axis of the dark 
matter. All of these facts can be understood in terms of 
the effects of angular momentum and dissipation during 
the merger. If there is a class of elliptical galaxies that 
were formed by gas-rich binary galaxy mergers, they can 
be expected to display these characteristics. 



Real galaxies in a ACDM universe are thought to have 
experienced many mergers over the course of their his- 
tory, and these multiple mergers can be expected to 
weaken the relationship between the shapes of galax- 
ies and their halos described here. The extent to which 
the effects of large-scale structure, such as mass accre- 
tion along filaments, tend to preserve the relationship 
between galaxies and their halos is an interesting and 
open question. 
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